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1. Introduction 

Experiments have raised the issue of the relative stability of various fractional quantum 
Hall states in graphene. The earliest results in suspended graphene in the two-terminal 
geometry [TJ gave indications of fractional quantum Hall effect (FQHE) at v = 1/3. 
The more recent experiments by Dean et al and Ghahari et al [2] in four-terminal Hall 
bar geometry for graphene on boron nitride hexagonal substrate have revealed a much 
richer structure. Specifically, they observe FQHE at filling factors \u\ = 2/3,4/3 and in 
the central (n = 0) Landau level, and at \v\ = 2 + 1/3,2 + 2/3,2 + 4/3,2 + 5/3,2 + 7/3 in 
the |n| = 1 Landau level. In addition, they also see evidence FQHE at \v\ = 1/3 in the 
n = LL, but none at \v\ =2 — 1/3. Hansel, Lafkioti, and V. Krstic have also reported 
4/3 FQHE in a four-terminal geometry [3|. More recently, Feldman, Krauss, Smet 
and Yacoby [4] have observed incompressible states in suspended flakes in the sequence 
v = p/(2p±l) for ah> < 4 in the \v\ < 1 range, but only at v = 2-2/3,2-2/5,2-4/7 
and 2 — 4/9 in the 1 < \v\ < 2 range. The difference of the < < 1 and 1 < < 2 
interval, and the absence of FQHE at 2 — 1/3 in particular, remains a puzzle. 

Our treatment below builds upon our understanding of the FQHE in GaAs systems. 
In GaAs systems, the FQHE is associated with the formation of composite fermions 
(CFs), bound states of electrons and an even number of quantized vortices j5]. The 
electronic LL splits into Landau-like levels of composite fermions, called A levels (ALs), 
and the integer quantum Hall effect of composite fermions produces the prominently 
observed FQHE at v = m/(2pm±l), where m is the number of completely filled A levels 
and 2p is the number of vortices bound to electrons. The CF theory has been extended 
to include spin [61 [7]: with spin, FQHE occurs at the same fractions as before, namely 
v = m/(2pm±l), but now, in general, FQHE states with different spin polarizations are 
available at each fraction, analogous to the situation at integer fillings where different 
spin polarizations would occur depending on the strength of the Zeeman splitting relative 
to the cyclotron energy. Experiments have measured the transitions between differently 
polarized states, the actual spin polarizations of the states, CF A level fan diagram, and 
the phase boundaries [8] , which are all found to be in good agreement with the CF theory. 
(For the integer quantum Hall effect, this physics is not relevant for GaAs, because the 
Zeeman energy is very small compared to the cyclotron energy, thus producing the state 
with the smallest spin polarization. However, for FQHE, the Zeeman energy can be of 
the order of the effective cyclotron energy of composite fermions, and therefore many 
different spin polarization occur.) We note that some of the states were written down 
in an early work of Halperin on multicomponent wave functions for the FQHE [9]. 

The FQHE problem in graphene differs from that in GaAs in two respects. First, in 
graphene, each electron has four-components, because of two spin projections and two 
valleys, producing an approximate SU(4) symmetry when the Zeeman energy and the 
valley splittings are negligible. Second, the linear dispersion leads to an interaction 
that is in general different from that in GaAs jTTJJ, [TTJ [12l [13l EH]. Our previous 
theoretical work on FQHE in graphene [TTJ [131 US] has shown that composite fermions 
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are formed in both n = and |n| = 1 LLs. The FQHE is thus still predicted to occur 
at v = mj (2pm ±1), where v = v — {An — 2) is the fractional part of the filling factor 
in LL n, but, in general, many more kinds of incompressible states are now available 
at each fraction, depending on whether the system is in the SU(2) limit or SU(4). The 
actual ground state is determined by a competition between the CF kinetic energy, 
which favors occupation of the lowest available A levels, the CF exchange energy, which 
favors occupation of A levels with the same SU(4) index, and the Zeeman and valley 
splittings [T3] . 

Some of the observations on FQHE in graphene are consistent with theoretical 
expectations: (i) The observed fractions have the form m / (2pm±l) as measured relative 
to a filled Landau level, (ii) The strength of FQHE in the |n| = 1 LL is comparable to 
that in the |n| = LL, to be contrasted with GaAs systems jTIJ [12], where FQHE in the 
second LL is much weaker than that in the lowest LL. (iii) As predicted in p3|[T5| [T6j. no 
FQHE has been observed at even denominator fractions in the \n\ = 1 LL. Some other 
observations are unexpected, however. Noteworthy is the absence of fractions such as 
2 — 1/3 and 2 — 2/5, while 2 — 2/3 and 2 — 4/3 are strong. Also puzzling is the appearance 
of FQHE at 1/3; one would have expected 2 — 1/3 to be stronger given that its reference 
state at v — 2 does not require any spontaneous symmetry breaking. In contrast, in 
the \n\ = 1 LL, 2 + 1/3 and 2 + 2/3 are both observed strongly. These fractions do 
not appear, at first sight, to be consistent with either SU(2) or SU(4) symmetry, thus 
posing a well defined puzzle whose resolution will shed important light on the nature of 
the FQHE in graphene. 

Our aim in this paper is to determine the stability of various states by evaluating 
the gaps using the CF theory. The gaps in the SU(2) limit were evaluated previously 
for states of the form v = m/(2pm + 1) for the lowest LL of GaAs (these results carry 
over to the n = Landau level of graphene [IT]) but not for the n = 1 Landau level 
of graphene. In addition, gaps for the reverse-flux attached states at v = p/(2pm — 1) 
have not been calculated for technical reasons. (Lowest Landau level projection, which 
is needed to obtain states appropriate for the high field limit, is prohibitively expensive 
for large systems for the latter class of states.) However, the reverse-flux attached states 
are relevant to the graphene problem because the observed 2/3 FQHE state is of that 
form [6j. We estimate below the gaps for the states relevant to the recent experiments. 
In addition, we obtain the collective mode dispersions for various graphene FQHE states 
in both the SU(2) and SU(4) limits. 

The issue of the 1/3 state in graphene has been addressed recently in exact 
diagonalization studies [18]. Another approach for constructing FQHE states in the 
SU(4) limit can be found in [19]. A Chern-Simons treatment of composite fermions in 
the SU(4) limit is given in [20] . 

Our paper is organized as follows. Following necessary definitions in section [2J 
we study the excitation structure of the SU(2) limit in section [3] and the SU(4) limit 
in section [H We discuss our results in light of the recent experiments j2j [3]. Our 
conclusions are presented in section [5j 
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2. Model 

2.1. Hamiltonian 

For electrons confined to a Landau level, the kinetic energy is quenched, and we need to 
consider the Hamiltonian that is diagonal in sublattice components (suppressed in our 
notation), 

Here e is the dielectric constant of the background material, &i and fj are Pauli matrices 
related to spin and Dirac valley indices, and Az and Ay are Zeeman and valley 
splittings. For Az = and Ay = 0, this Hamitonian possesses an SU(4) symmetry, with 
<7 Z C§)1, l®f z , o~ z ®f z (or any appropriate linear combinations thereof) defining the weights 
in SU(4) multiplets. Composite fermions acquire additional labels (quantum numbers) 
cr 1 ), . . . , o^ 4 -*, which come from a basis chosen in the fundamental representation of 
SU(4), and can be chosen conveniently to be t1\ T4j It an d 11 where the first arrow 
refers to valley and the second to spin. Let (777.1,777.2,777.3,777.4) denote the incompressible 
ground state in which 777.1, • • • > m 4 A levels with label . . . , are fully occupied. 
Which of these states is the ground state depends [13] on the values of A z and Ay. A 
coupling between the spin and valley degrees of freedom and LL mixing are neglected 
in what follows. 

We neglect the anisotropy of the interaction due to underlying sublattice structure 
of the basis states, because such effects are small, being proportional to the ratio of the 



lattice constant a to the magnetic length £ = ^jTic/eB. The study of the residual lattice 
effects is beyond the scope of this paper. 

Another important assumption of our model is the neglect of disorder. This is a 
standard assumption in most quantitative studies of FQHE, because it is not known how 
to include the effect of disorder in a reliable manner. It is possible that disorder might 
be affecting different FQHE states differently, thus affecting their relative stability. In 
the absence of a reliable quantitative treatment of disorder, while comparing our results 
with experiment we will make the simplest assumption that while the disorder causes 
a significant reduction of the excitation gap, it does not change the order of stability of 
various FQHE states. Strictly speaking, our results below ought to be considered as a 
prediction for future experiments in the limit of low disorder. 

The kinetic term still determines the Coulomb matrix elements via the single- 
particle orbitals when the Landau level is fixed. For the n = LL the interelectron 
interaction is same as in GaAs, because in this LL the single-particle states are identical 
for carriers with linear and quadratic dispersions. For n = ±1, however, the two- 
component character of the single-particle states of massless Dirac fermions produces 
different interaction pseudopotentials than in GaAs. Following standard practice, 
[TT1 [TBI [15] we use a lowest LL basis with an effective interaction given in [T3j , which 
reproduces the \n\ = 1 graphene pseudopotentials. We employ the standard spherical 
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geometry for our calculations [21], in which electrons move on the surface of a sphere 
and a radial magnetic field is produced by a magnetic monopole of strength Q (integer 
or half-integer) at the center, producing a magnetic flux 2Q(hc/e) through the surface. 

2.2. CF wave functions 

The general form of the CF wave functions in the spherical geometry is given by 

^ = V LLL & m J] (ujV k - v jUk ) 2p , (2) 
j<k 

Here, u = cos (9/2)e~ i ^ 2 and v = sin(#/2)e^/ 2 . tf„ is the wave function for interacting 
electrons at v and $ m is the wave function for noninteracting electrons at filling factor 
m, with v = m/(2pm ± 1). In the spherical geometry, \l/ m is constructed at monopole 
strength q, producing \1/ at monopole strength Q = q + p(N — 1). "Plll projects the 
state into the lowest LL [22] . The orbital part of the wave function has the same form 
as above even for the SU(2) and SU(4) generalizations, which we discuss below. 

2.3. Effective interaction 

The electron-electron interaction is conveniently described in terms of pseudopotentials 
[2~T] V m , which give the energy of two electrons in relative angular momentum m. As 
the orbitals in the lowest LL in graphene are identical to those in the conventional 
two-dimensional electron gas, the pseudopotentials are identical too. The problem of 
interacting electrons in the n = 1 LL of graphene can be mapped into a problem of 
interacting electrons in the lowest LL with the effective pseudopotentials [12] 

^-/StK £ '(t)^(t))>^ < 3 > 

where L n are Laguerre polynomials. 

The Monte Carlo evaluation of the energy of variational wave functions, on the 
other hand, requires a real-space interaction. To simplify calculation, we construct the 
wave functions in the lowest LL, and use an effective real-space potential corresponding 
to the effective pseudopotentials V m in the lowest LL basis. Following a well-tested 
procedure [231 H]> we use the form 

w 1 M 

V eS (r) = - + J2c i r i e- r , (4) 

and fit the coefficients q to reproduce the first M+l pseudopotentials Vm 'of equation 
([3]). The coefficients Cj obtained for M = 6 are given in [13]. 

3. SU(2) analysis 

To weigh the relative importance of the relevant parameters, it is appropriate to consider 
the ratio of the Zeeman splitting to some relevant energy scale in the problem. We 
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consider « = A z /(0.02e 2 /<) = Q.Q5gey/B[T\, where the energy in the denominator is 
the Fermi energy of the single component CF Fermi sea and g is the Lande factor for the 
background semiconductor material [7J [8] . (The quantity £ = ^Tic/eB is the magnetic 
length. We could alternatively have chosen the gap of a prominent FQHE state.) The 
FQHE states in GaAs are typically fully spin polarized for B > 5 T, i.e., for ratio 
k >0.5 [TJ IB]. For graphene, this ratio is 1.8 (0.6) at B = 36 T (4 T), suggesting that 
the spin degree of freedom may be frozen in the experiments of [2], which are performed 
at rather large B. The valley symmetry breaking in graphene is on the order of a /I, 
where a is the lattice spacing [T4J. At B = 4 and 36 T, this ratio is a/l pa 0.02 and 
0.06, respectively. While this issue requires further investigation, it cannot be ruled 
out that the spin is fully polarized but the valley symmetry may not be broken [24J. 
We therefore first consider the situation where the spin degree of freedom is frozen but 
valley symmetry is intact. (The analysis of the SU(4) limit is given in the next section.) 

In the SU(2) limit, the orbital part of the wave functions at the maximal weight 
states [25] of an SU(2) multiplet at monopole strength Q are given by equation (T5]), with 




with mi + m 2 = m and the superscript denotes the electron species. In other words, the 
state consists of m\ and m 2 Landau bands of up and down pseudospins, respectively 
occupied. (Here, pseudospin refers to the valley SU(2) symmetry.) One can show that 
this wave function satisfies Fock's cyclic condition |25j, i.e., it is annihilated by any 
attempt at further antisymmetrization; this ensures that it is an eigenstate of the total 
pseudospin operator. The energies of variously valley-polarized states can be evaluated 
using standard methods, from which the phase diagram as a function of the valley 
Zeeman energy can be constructed [7J. The valley degree of freedom behaves very 
much like the spin degree of freedom. In fact, the predictions of the CF theory for two 
component systems have been verified semiquantitatively [26] for another system with 
two valleys, namely AlAs quantum wells, where the valley splitting (analogous to the 
Zeeman energy) can be controlled by application of in-plane symmetry breaking strain. 

With spin frozen, 1/3 is valley-polarized, whereas 2/3 and 2/5 are valley singlets. 
(These states are related by particle-hole symmetry to 5/3, 8/5, and 4/3.) That 9 = 2/5 
and 2/3 are singlets in the n = Landau level follows from the analogous result 
established in conventional two-dimensional electron systems such as GaAs |6] . 

For the v = 2/5 in the |n| = 1 LL of graphene the energies of the 2/5 FQHE 
states in the thermodynamic limits can be obtained using the CF wave functions, 
and are given by -0.4450(2)e 2 /e£ for the singlet and -0.4401(2)e 2 /e£ for the polarized 
ground state energies, as shown in Fig. HJa). The polarization energy, i.e., the energy 
difference between the polarized and the SU(2) singlet states, has been estimated to be 
0.00497(9)e 2 /e£ in the thermodynamic limit as shown in Fig. Hfb). This confirms that, 
to the extent the wave functions are accurate, the energy ordering of states predicted by 
the CF theory remains same in the |n| = 1 LL as in the n = LL. We note that Shibata 
and Nomura [27] conclude, based on a density matrix renormalization group calculation, 




m 




(5) 
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Figure 1. (a) The extrapolation of the energy of the fully polarized and the SU(2) 
singlet state at v = 2/5 in the \n\ = 1 Landau level of graphene to the thermodynamic 
limit. N is the number of particles. Only the leftmost four points, corresponding to 
N > 48, are included in the fit because smaller systems show finite size effects, (b) 
The polarization energy at v = 2/5. The scaling is linear for N > 32. 



that the ground state here is polarized even in the absence of a symmetry breaking field; 
this conclusion, different from ours, is based on a study of smaller systems and without 
extrapolation to the thermodynamic limit. 

We assume in what follows that the valley splitting is small but large enough 
to suppress the CF valley- skyrmions (CF skyrmions are rapidly suppressed by the 
Zeeman/valley splitting [28]), and the relevant excitation is a particle- hole pair of 
composite fermions. If composite fermions were completely noninteracting, the 1/3 
gap would be the smallest, as it does not involve CF cyclotron energy. The composite 
fermions, however, do feel a residual interaction, which necessitates a quantitative 
evaluation of the gaps, which can be accomplished using the standard methods of the 
CF theory [22| [29]. The energy gaps of the lowest charged excitations shown in Tabled] 
have several unexpected features [30]: The 2/3 and 4/3 states in the n = 1 LL are the 
strongest, followed by 2/3 and 4/3 in the n = LL and 1/3 and 5/3 in the n — 1 LL. 
Next most stable are 1/3 and 5/3 in the n = LL. These calculations are consistent 
with the observations that 2/3 and 4/3 are the strongest in the n = LL (producing 
FQHE at 2 — 2/3 and 2 — 4/3), both 1/3 and 2/3 are strong in the n = 1 LL (resulting 
in FQHE at 2 + 1/3, 2 + 2/3, 4 ± 1/3, 4 ± 2/3), and 2/5 and 8/5 are not yet observed. 

While the observed order of stability is consistent with theory, the actual values of 
gaps are not. The largest observed gap [2]of~16Katz/ = 4/3 and B = 35 T is much 
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v 



n = 



n\ = 1 




0.073(1) 
0.11(4) 
0.0390(6) 



0.095(2) 
0.21(14) 
0.0475(9) 



Table 1. The smallest gaps for the creation of a far separated particle-hole pair of 
composite fermions. The 1/3 ground state is fully spin and valley-polarized, whereas 
the 2/3 and 2/5 states are spin polarized but valley singlets. Neglecting LL mixing, 
the gaps at v and 2 — v are equal. All energies are quoted in units of e 2 /e£. The 2/3 
gaps come from a combination of exact diagonalization and CF theory; the rest comes 
from CF theory only. 



smaller than the theoretical estimate ~ 0.1e 2 /e£ ~ 110 K (for e = 3). The significant 
discrepancy is no doubt in large part due to disorder, and suggests that improvements 
in sample quality should reveal much further FQHE structure. 

4. SU(4) analysis 

We next turn to the SU(4) limit, achieved when the valley and the spin states are all 
approximately degenerate and may be applicable to experiments at very low B. In this 
limit, if composite fermions were completely non-interacting, only FQHE states of the 
form m/(2pm±l) with m = 4, 8, • • • would be observable, but as seen above, the residual 
interaction between composite fermions results in a spontaneous symmetry breaking, 
producing FQHE at all of these fractions [13] . We now show that new Goldstone modes 
become available, which are an SU(4) generalization of the spin waves of composite 
fermions considered previously in the context of GaAs FQHE [31]. (Similar mechanism 
has been considered for QHE at v ^ An — 2 integer fillings in graphene [ID] and in bilayer 
systems [32].) As a result of the presence of these modes, fractions such as 2/5 and 3/7 
are expected to be much weaker than their SU(2) counterparts. 

In the SU(4) generalization [13] of CF theory [191 ED], the orbital part of the ground 
state wave function at v = m/ (2pm ± 1) for correlated electrons in the maximal weight 
states [25] of an SU(4) multiplet at monopole strength Q is still given by equation El 
but with 



at monopole strength q = Q — p(N — 1), with iV = J2 S N S and m = J2 s m s- With 
mi > m2 > m 3 > m i) it is easy to see that \1/ satisfies Fock's cyclic condition 
[25], i.e., it is annihilated by any attempt to antisymmetrize an electron / of type u 
(min M < / < max u ) with respect to the electrons of type t < u, 





(7) 
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where (k, I) permutes indices k and /. For the CF-particle hole excitation, one must in 
general choose $ to be an appropriate linear combination of Slater determinant states 
to ensure the Fock condition. In each SU(4) multiplet, represented by a Young tableau 
(YT) [Mi, M 2 , M 3 ], where M t is the length of the z-th row (M« > M i+1 and the last M^s 
are not shown if zero), the maximal weight states can be constructed with the help of 
elementary group theory, and the eigenstates of the total orbital angular momentum 
L are obtained by standard Clebsch-Gordan expansion. We show the appropriate 
combinations pictorially for all single exciton excitations at 2/5, 3/7 and 4/9, along 
with their YT representations, in the subsequent figures. The interaction energies of 
the explicit wave functions are evaluated by Monte Carlo integration. For CF particle- 
hole excitation we identify their orbital angular momentum L with their wave vector k 
according to k — L/£, which yields the dispersion for the excitations, with the k — > oo 
limit giving the energy of a far separated CF particle-hole pair. All states in a given 
multiplet are degenerate for A z = = Ay; for nonzero A z and Ay the energy splittings 
are straightforwardly determined. 

The CF ground state at v — 1/3 has a single fully occupied n = AL. As excitation 
to an unoccupied n = AL is already possible in SU(2) systems, we expect no new 
physics at this fraction. 

The v = 2/5 state has two n = ALs occupied in the ground state; assuming 
for defmiteness that Az > Ay, the occupied ALs are tt an d ti- It nas two kinds of 
excitations depicted in Fig. [2j (i) The highly degenerate SU(4) spin wave, in which 
a CF is excited "sideways" to an unoccupied n = AL, as shown in Fig. [21(a)- This 
mode, represented by the YT [N/2, N/2 — 1, 1], is available only because of the SU(4) 
symmetry, (ii) Excitations within the SU(2) subspace of the ground state, in which 
a CF is raised "up" to the second AL either flipping or preserving its valley index, 
represented, respectively, by YTs [A^/2 + I, N/2 — 1] and [N/2, N/2]. The four cases 
give rise to a valley triplet (Fig. Mjo)) and a valley singlet (Fig. EJ^c)), which are both 
gapped at all wavelengths. The triplet is split by the Ay and destabilizes the FQHE state 
when it reaches the "valley-roton" gap [31] of Agfg 1 = 0.026(5)e 2 /el The energy of the 
"singlet-roton" gap in mode (c) is unaffected by valley or Zeeman splittings. Transitions 
that raise a CF to the first excited level of an otherwise empty band are ignored, as 
these involve a loss of exchange energy and a kinetic energy cost at the same time and 
are likely to represent higher excitations. 

Fig. [2] shows that the smallest gap to creating a far-separated pair of charged 
excitations, which is the one relevant for transport experiments, corresponds to a CF 
exciton in mode (a). Because this energy is approximately half the gap in GaAs for 
either fully spin polarized or spin singlet 2/5 state [31], the new Goldstone mode thus 
results in a substantial weakening of the 2/5 state in the SU(4) limit. From the SU(4) 
spin wave dispersion at small k we estimate the energy of the SU(4) skyrmion, using 
the methods of [33J, to be 0.030(6) and 0.037(9) in the n = and \n\ = 1 Landau levels; 
this has a higher energy than the other excitations, so will not be relevant to either 
transport or spin polarization. 
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Figure 2. Excitations and gaps of the v = 2/5 fractional quantum Hall effect 
in the \n\ =0,1 graphene Landau levels. The A level structure for the maximum 
weight state in each mode is indicated in the diagrams at the top. Mode (a), which 
has no analog in GaAs, determines the smallest charge gap. The modes (b) and (c) 
correspond to singlet and triplet excitations in an SU(2) sector embedded in a SU(4) 
multiplet. The transport and roton gaps shown near the largest momenta and the 
minimum, respectively, correspond to the thermodynamic limit. The rotons in mode 
(b) destabilize the state when Az or Ay becomes greater than the roton energy. 



At v = 3/7 there exist "sideways" excitations (Fig. E(a)) as well as those that 
stay within the SU(3) sector of the ground state (Fig. [3] (b) and (c)). The nine possible 
excitons within the SU(3) sector can be reduced to a [2, 1] octet of SU(3) (see illustration 
above Fig. Mb)) and an SU(3) singlet. The latter SU(3) multiplets are embedded in 
[N/3 + 1,2, 1] and [N/3] SU(4) multiplets, respectively, and are all gapped. The SU(4) 
spin wave mode of Fig. [3(a) is not gapped in the k — > limit and has the smallest energy 
at all wave vectors. The rather small gap to creating a far separated CF particle-hole 
pair in mode (a) makes the 3/7 state very delicate in the SU(4) limit. The stiffness of 
the SU(4) spin wave in (a) could not be determined with sufficient accuracy, precluding 
an estimate of the SU(4) skyrmion energy at 3/7. 

Weak symmetry breaking fields A z and Ay split the multiplet and select a state of 
reduced symmetry out of the SU(2) xSU(2) decomposition [M] of the SU(4) multiplet for 
both the ground state and the exciton. To illustrate we consider 9 = 2/5 with Az > Ay. 
The spin and pseudospin (valley) quantum numbers are (S, P) = {N/2, 0) for the ground 
state and (N/2 - 1, 1 or 0), (N/2,1) and (JV/2,0), respectively, for modes (a) to (c). 
Thus the fourfold degenerate mode (a) is raised by A z and split by —Ay, 0, 0, Ay; mode 
(b) is split similarly but it is not shifted; mode (c) is neither split nor shifted. For each 
mode as well as the ground state there are other copies displaced to higher energies in 
multiples of Az, which in turn must merge into the continuum of other excitations; the 
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Figure 3. Excitations and gaps [29] of the v = 3/7 fractional quantum Hall effect 
in the |n| =0,1 graphene Landau levels. The spin-wave-like mode (a) produces the 
smallest charge gap. Modes (b) and (c) produce an SU(4) roton; the roton in (b) causes 
an instability with increasing Az or Ay. Roton and transport gaps are indicated. 2x 
denotes weight multiplicity. 

structure of these higher copies follows from the SU(4) to SU(2) xSU(2) decomposition 

We have also studied v = 4/9, where there is no spin wave but the sixteen possible 
excitations into the second AL decouple into an SU(4) singlet (b) and the fifteen-fold 
degenerate mode (a) labeled by the YT [2, 1, 1], see Fig. HJ The large k limits of lowest 
energy are 0.021(1) and 0.026(2) for the n = and n = 1 LLs. As a function of Az 
or Ay there are several possible transitions into CF states with other spin and valley 
quantum numbers [13]. The roton and transport gaps in mode (b) could not be obtained 
with sufficient accuracy in the theormodynamic limit; but in any finite systems they are 
not smaller than in mode (a). It is uncertain if mode (b) truly develops a roton gap. 
Rotons in mode (a) destabilize the state at sufficiently high Zeeman energy. 

5. Conclusion 

To conclude, we have argued that the spin degree of freedom may be frozen at high fields 
in graphene, and that many features of the FQHE in graphene at high magnetic field are 
consistent with SU(2) CF theory. In particular, quantitative calculations are consistent 
with the observed order of stability of states in n = and |n| = 1 LLs, although the 
observed gaps are much smaller than the theoretical values. We have further shown 
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Figure 4. Excitations of the SU(4) singlet v = 4/9 FQHE in the |n| = 0, 1 Landau 
levels. The weight structure of the 15-fold degenerate mode (a) is shown. Notice the 
roton minimum in mode (a), which is destabilizing if Az + Ay ~ A^~ g ro * on . 3x 
denotes weight multiplicity. Roton and transport gaps in mode (a) are indicated. 



that in the SU(4) limit, the FQHE states at v = 2/5 and 3/7 have a greatly reduced 
charge gap due to new SU(4) spin wave mode, which implies more stringent conditions 
for their observation. 
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